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1 Introduction 

Models and event generators for hadronic collisions belong to the basic building 
blocks of hadron cascade calculations to simulate radiation problems around high 
energy accelerators and to simulate the Cosmic Ray cascade in the atmosphere, 
within the earth or within space stations. At the same time these models often 
have to describe also the new physics to be investigated at a given accelerator, 
examples are the present study of photon-photon collisions at the LEP collider 
or the study of central collisions of heavy ions at the CERN-SPS and at the 
RHIC and LHC colliders. 

There are at least two main applications of such models to study radiation 
problems around accelerators and colliders: 

(i) At colliders like the LHC one of the main radiation sources are the proton- 
proton or nucleus-nucleus collisions inside the detectors. Present and future 
colliders are also to investigate e + - e - , /i + - fj,~ and photon-photon collisions. 
Therefore, also models for hadronic collisions of leptons and photons become 
more and more important for the study of radiation problems. 

(ii) The second application is for the transport of the produced particles in 
the elements of the detectors and accelerators and in the shielding material. For 
this we have to understand hadron-nucleus, nucleus-nucleus, lcpton nucleus, 
photon-nucleus and even neutrino-nucleus collisions at all energies. 

Here we are concerned with hadronic models or event generators for hadronic 
collisions. These models are inserted into hadron cascade codes to study the 
radiation transport. We will not discuss hadron cascade codes here, but some 
of them are mentioned when discussing the application of the hadronic models. 
There are also cases, where a dedicated hadronic model is formulated within the 
hadron cascade code. 

The models for hadronic collisions have to provide the produced hadrons 
preferably in the form of Monte Carlo events as well as the inelastic cross sections 
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of the collisions considered. The event generators of interest for our applications 
have to provide minimum bias events, not only events for selected (hard) colli- 
sions, which might be the main object of the physics studied at a given collider. 
We need detailed models for all collisions, not only for interesting events, where 
for instance SUSY particles or Higgs bosons might be produced. In recent years 
new features of hadronic collisions have been found experimentally which were 
not predicted by any of the models. Therefore, the models have to be contin- 
uously updated and revised in order to reproduce all known properties of the 
collisions. Some of these features are: (i)Hard diffraction and events with 
large rapidity gaps between jets J|,(| . (ii) Enhanced stopping of leading baryons in 
nuclear collisions ^,|) , one way to incorporate this into the multi-string models 
is by considering new string configurations (diquark-breaking diagrams) 0,0]. 

2 Models for hadronic interactions 

Unfortunately, it is not possible to discuss here all existing Monte Carlo models 
for hadron production. Let me at least mention some important models, which 
I will not discuss: 

(i) There are JETSET §, PYTHIA fo), HERWIG gl) and IS A JET @, 
which are mainly used to sample selected classes of events. These are certainly 
the event generators used most extensively by collider experiments. Codes not 
discussed, mainly for reasons of space, are also HIJING (H and FRITIOF @, 
both are codes for nuclear collisions. 

(ii) There are dedicated event generators inside hadronic codes like FLUKA, 
HERMES, MARS, MCNPX, NMTC or GEANT. Usually these are not available 
as separate codes. The event generators inside all of these codes will be covered by 
contributions to this Meeting (examples: FLUKA: |]| GEANT: pi). Therefore, 
there is no need for me to discuss them. 

(iii) PHOJET is a minimum bias Dual Parton Model event generator 
with a very broad applicability for hadrons, photons and electrons as projec- 
tiles and targets. PHOJET is contained in the new DPMJET-III event gener- 
ator, which adds nuclei as targets and projectiles to the above PHOJET list. 
DPMJET-III will for the first time be described at this meeting JT§ . 

In this talk I will only discuss hadronic models which work above a mimimum 
lab-energy of about 5 GeV per nucleon and with a maximum lab-energy of about 
10 11 GeV or (in the units used in Cosmic Ray applications 10 21 eV), this means 
codes which are also available for studying the Cosmic Ray cascade up to the 
highest energy. 

The extension of models for multiparticle production in hadron- hadron, 
hadron-nucleus and nucleus-nucleus collisions to be used for the simulation of 
the Cosmic Ray cascade up to Ei a \, = 10 21 eV (corresponding for p-p collisions 
to y/s — 2000 TeV) is needed to prepare for the Auger experiment pof as well 
as for a reliable interpretation of present experiments like AGASA Fly's 
Eye |Q and High Resolution Fly's Eye [p3|j2^| , which present data in the EeV 
energy region. 
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Nearly all event generators for hadronic or nuclear collisions in the considered 
energy range are constructed out of the following ingredients: 

(i) In hadron-hadron collisions or in elementary nucleon-nucleon collisions of 
nuclear collisions multiple soft color singlet chains with quarks, diquarks and 
antiquarks at their ends are formed according to the rules of Gribov-Regge 
theory together with topological arguments. This are mainly diquark-quark and 
quark-antiquark chains. The Gribov-Regge theory also gives all the hadron- 
hadron cross sections. 

(ii) In collisions involving nuclei one starts with the construction of multiple 
interactions according to the Gribov-Glauber theory. In simple models each 
collision leads just to the formation of one pair of colorless chains. In evolved 
models there is the full multi-chain system according to (i) and (iii) in each 
elementary collision. The Gribov-Glauber theory gives not only the distribution 
in the numbers of elementary collisions, it gives also the nuclear cross-sections 
in terms of the elementary (input) nucleon-nucleon cross-sections. There are 
also simplified models, which do not treat the full Gribov-Glauber theory for 
nucleus-nucleus collisions. In these models (using the superposition model) the 
nucleus-nucleus collision is replaced by a couple of hadron-nucleus collisions. 

(iii) At sufficiently large energy besides the soft chains according to (i) and 
(ii) we have a system of multiple minijets. At high collision energies hard and 
semihard parton-parton collisions occur and become more and more important 
with rising collision energy. The produced jets and minijets become important 
features of most models, they determine the transverse momentum or trans- 
verse energy properties of the models at high energy. The input cross section 
for semihard multiparticle production (or minijet production) Oh is calculated 
applying the QCD improved parton model, the details are given (for instance as 
implemented in DPMJET) in Ref . |^M|7pJ] . 
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fi(x,Q 2 ) are the structure functions of partons with the flavor i and scale Q 2 
and the sum i,j,k,l runs over all possible flavors. Now after the HERA mea- 
surements of structure functions at small Bjorken Xbj it is very essential to use 
modern structure functions, which agree to the HERA data. Since the HERA 
measurements, the structure functions are known to behave at small x like l/x a 
with a between 1.35 and 1.5. To remain in the region where perturbation theory 
is valid, a low p± cut-off, p± thr , is used for the minijet component. The multi- 
plicity distribution for the minijets is obtained by an unitarization procedure. 

(iv) After forming all partonic chains according to (i) to (iii) the next step 
is initial state and final state QCD evolution of all hard chains, after this all 
chains are hadronized using the Lund codes JETSET or PYTHIA or in some 
codes private dedicated codes for string hadronization. 

(v) There are special cases of hadron production in hadron-hadron, hadron- 
nucleus and nucleus-nucleus collisions, these are diffractive interactions (or in- 
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teractions with rapidity gaps between the produced particles) . Diffractive hadron 
production is usually treated in special routines. Since the discovery of hard pro- 
cesses also in diffractive interactions at the CERN-SPS collider, HERA and at 
the TEVATRON (for instance jet-gap-jet events) this has become more com- 
plicated than before, only very few codes (for instance PHOJET and with this 
DPMJET-III) include at present hard diffraction. However, there exist dedicated 
codes for the simulation of hard diffraction. 

(vi) After steps (i) to (v) we have events consisting out of produced hadrons 
(from string fragmentation) and spectators (these are nucleons from the original 
projectile and target nuclei, which have not taken part in the interaction). In 
the next step the codes treat secondary interactions. This is in simple cases 
just a formation zone intranuclear cascade, where the produced hadrons interact 
with the spectators, a more complete treatment is the full secondary hadron 
cascade, where the produced hadrons interact with the spectators as well as with 
other produced co-moving hadrons. There are also codes, where the secondary 
interactions occur already at the level of the chains before hadronization or at 
the level of chain-end partons before the formation of the color neutral chains. 

(vii) Now we have the produced hadrons and the remains of the projectile and 
target nucleus (in the form of spectator nucleons with no interactions or only in- 
teractions which were not able to knock them out of the nuclear potential). These 
remains are considered as excited residual nuclei, which in a following evapora- 
tion step deexcite by the emission of evaporation particles (protons, neutrons 
and nuclear fragments). Finally the excited nuclei, which no longer are able to 
emit evaporation particles emit deexcitation photons to form finally a stable (or 
radioactive) residual nucleus. 

Let us summarize important aspects of some event generators used often in 
calculating the Cosmic Ray cascade: 

The DPMJET-II.5 event generator based on the two-component Dual Par- 
ton Model (DPM) was described in detail ^p^QHHg] . DPMJET-II.5 de- 
scribes well minimum bias hadron and hadron jet production up to present 
collider energies. DPMJET is used for the simulation of the Cosmic Ray cascade 
within the HEMAS-DPM code Q used mainly for the MACRO experiment 
[B5 . DPMJET-II.5 will soon be available within the CORSIKA cosmic ray cas- 
cade codep6[. 

The SIBYLL-1.7 model |37f is a minijet model and has been reported to be 
applicable up to Ei a b = 10 20 eV. However, the EHLQ (38) parton structure func- 
tions used for the calculation of the minijet component might , after the HERA 
experiments, no longer be adequate. It is known, that a significant updating of 
SIBYLL is on the way |39), but the new code is not yet available. SIBYLL is 
a popular model for simulating the Cosmic Ray cascade in the USA. SIBYLL 
is available within the CORSIKA cosmic ray cascade code|36| and within the 
AIRES go| and MOCCA gl],|| cosmic ray cascade codes . 

VENUS f43[| , a popular model applied originally mainly for describing heavy 
ion experiments. VENUS is applicable up to E^ — 5xI0 16 eV. It has been 
reported [Q, that the introduction of minijets into VENUS has been planned, 
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this will allow to apply VENUS up to higher energies. A new version of VENUS 
called neVfis[^5| was not yet to be available for distribution at the time, when 
my talk was prepared, it will be available in the next release of the CORSIKA 
cosmic ray cascade code. 

QGSJET Q is the most popular Russian event generator used for Cosmic 
Ray simulations. It is based on the Quark Gluon String (QGS) model, this model 
is largely equivalent to the DPM. QGSJET also contains a minijet component 
and is reported to be applicable up to Ei a b = 10 20 eV. QGSJET is available 
within the CORSIKA cosmic ray cascade code|3(| and within the AIRES j4(j| 
cosmic ray cascade code. 

In Table 1 we present some characteristics of the models. The Gribov-Regge 
theory is applied by three of the models. The pomeron intercept for SIBYLL 
is equal to one, SIBYLL is a minijet model using a critical pomeron, with only 
one soft chain pair, all the rise of the cross section results from the minijets. In 
the models with pomeron intercept bigger than one, we have also multiple soft 
chain pairs, already the soft pomeron leads to some rise of the cross sections with 
energy. Minijets are used in three of the models, it is believed, that minijets are 
necessary to reach the highest energies. All models contain diffractive events. 
Secondary interactions between all produced hadrons and spectators exist only 
in VENUS, DPMJET has only a formation zone intranuclear cascade (FZIC) 
between the produced hadrons and the spectators. Only three of the models 
sample properly nucleus-nucleus collisions, SIBYLL replaces this by the semi- 
superposition model j47|. The residual projectile (and target) nuclei are only 
given by two of the models. 

Table 1. Characteristics of some popular models for hadron production in 
Cosmic Ray cascades. 
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Each model has to determine its free parameters. This can be done by a 
global fit to all available data of total, elastic, inelastic, and single diffractive 
cross sections in the energy range from ISR to collider experiments as well as 
to the data on the elastic slopes in this energy range. Since there are some 
differences in the hard parton distribution functions at small x values resulting 
in different hard input cross sections we have to perform separate fits for each set 
of parton distribution functions. After this stage each model predicts the cross 
sections also outside the energy range, where data are available. 



3 Code comparisons 

There is certainly the need to understand the systematic errors of the hadronic 
models and event generators, (i) This is more easy in all the region, where 
experimental data exist. Here it is useful to compare the models to as many 
experimental data as possible, ft is a task of the community to decide about 
the most suitable benchmark experiments. We can however assume, that most 
published models agree to most of the available data. 

Systematic errors of models are more difficult to estimate in regions, where 
no experimental data are available. This might be collisions at energies beyond 
the energies of existing accelerators or in regions of the phase space, where 
no experimental data are available. A very prominent example for the latter 
is the fragmentation region at large Feynman-x, which was never studied in 
hadron production experiments at the highest energies available at the CERN- 
SPS collider and at the TEVATRON collider, but just this fragmentation region 
is of upmost importance for the Cosmic Ray cascade. The systematic comparison 
of Monte Carlo models in such regions is very useful to understand how reliable 
the extrapolations of the models into such regions might be. 

The Karlsruhe code comparison |Q was a first extensive comparison of 
Monte Carlo models, which are available as event generators within the COR- 



SIKA Cosmic Ray cascade code |36|j49| . The energies, mainly in the 10 14 and 
10 15 eV range and the distributions chosen in this comparison were motivated by 
the interest of the KASKADE|5(| experiment in Karlsruhe. With one exception, 
I will not review this code comparison here. A code comparison extending up to 
higher energies was performed within the AIRES jlO) code system | f3l| . 

Further code comparisons should be strongly encouraged. Some extensions 
to higher energy of such a code comparison were presented in talks by Engel 

[mi- 

Here I will present some code comparisons up to the highest energies, for 
which the models are able to run. The first idea was just to use the event gen- 
erators implemented in the last distributed version CORSIKA-5.61 |49| . But 
there was only the older version DPMJET-II.4 implemented. So I use for DP- 
MJET my stand-alone version of DPMJET-II.5 ||||,@]. With SIBYLL-1.6 as 
implemented in CORSIKA-5.61 I got funny results at high energy (above 10 18 
eV) (I understand, that these problems with SIBYLL-1.6 are now corrected in 
CORSIKA.). Using the stand-alone code SIBYLL-1.7, the last SIBYLL ver- 
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sion distributed by the authors, these anomalies disappeared, therefore I use for 
SIBYLL only the stand-alone version mentioned. With QGSJET in CORSIKA, 
I did not find any problems, so I use it in the form as implemented in CORSIKA- 
5.61. The VENUS code as implemented in CORSIKA-5.61 runs only up to 10 16 
eV. With this it is not possible to go beyond the Karlsruhe code comparison 




In all the following plots I present and compare the result of the three event 
generators in the implementations as mentioned. The first comparisons are for 
p-N (proton-Nitrogen) collisions, p-N collisions are a good approximation to 
p-Air collisions which are important for the Cosmic Ray cascade. 
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Fig. 1. (a) Average transverse momenta for charged hadron production in p-N 
(proton-Nitrogen) collisions as function of the (nucleon-nucleon) cms energy yfs. (b) 
Average transverse momenta for proton production in p-N (proton-Nitrogen) collisions 
as function of the (nucleon-nucleon) cms energy y/s. 



In Fig.|l|.a we present average transverse momenta of charged hadrons as ob- 
tained from DPMJET-II.5, QGSJET and SIBYLL for p-N collisions as function 
of the cms energy y/s. At energies where data (in p-p or p-p collisions) exist all 
models agree rather well with each other and with the data. However, we find 
completely different extrapolations to higher energies. We should note, all three 
models have a minijet component. But it seems, that in spite of the minijets the 
average transverse momentum in QGSJET becomes constant at high energies, 
while it continues to rise in DPMJET. We can conclude, there are very big dif- 
ferences in implementing the minijet components in the models. There are even 
bigger differences in the models, if we consider the average transverse momenta 
of secondary protons, see Fig.[j].b. DPMJET and QGSJET diverge like in Fig.|l|.a 
and the prediction of SIBYLL-1.7 is even at lower energies, where we find agree- 
ment between all three models in Fig.[l|.a, significantly outside the other models. 

Fig.|].a presents the rise of the total charged multiplicity and FigJ|.b presents 
the rise of the secondary proton multiplicity with the cms energy y/s according to 
DPMJET, QGSJET and SIBYLL. We find again, at low energies, where data are 
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Fig. 2. (a) The average multiplicity of charged hadron production in p-N (proton- 
Nitrogen) collisions as function of the (nucleon-nucleon) cms energy ,/s. (b) The av- 
erage multiplicity of secondary proton production in p-N (proton-Nitrogen) collisions 
as function of the (nucleon-nucleon) cms energy ^fs. 



available, the models agree rather well. DPMJET and SIBYLL agree in all the 
energy range shown. However, QGSJET above the energy of the TEVATRON 
extrapolates to higher energies in a completely different way. 

In Fig.||.a and ||.b we present for p-N collisions the energy fractions K of 
charged pion production and for net-baryon B — B production (Since newly 
produced baryons and anti-baryons are expected to carry the same energy frac- 
tions, this corresponds to the energy carried by the baryons which were present 
before the collision, which have been stopped and transformed by the collision.). 
In the case of K B _ B we find all three models to behave widely different, in the 
case of K n +^- we find DPMJET and QGSJET to agree largely, but SIBYLL 
behaves quite differently. Here I should mention, that new baryon stopping di- 
agrams [|4]|5f| have so far only been implemented in DPMJET-II.5. This en- 
hanced stopping is responsible for K B _ B being smaller and i7r - being bigger 
in DPMJET-II.5 than in the other two models. 

The cosmic ray spectrum-weighted moments [ |56| in p-A collisions are defined 
as moments of the F(xi a b) = xi a bdN/dxi a b ■ 

fT A = I {xia b T- l F?- A {x lab )dx lab (2) 
Jo 

Here —7 ~ -1.7 is the power of the integral cosmic ray energy spectrum and A 
represents both the target nucleus name and its mass number. 

In FigJ^.a and ^.b we present the spectrum weighted moments for pion and 
Kaon production in p-N collisions as function of the cms energy -^/s per nu- 
cleoli. In contrast to Fig. || we find this time DPMJET and SIBYLL to agree 
approximately, but QGSJET behaves completely different. 

In Figs. |[a and ||.b we present the extrapolations to high energy Ei ab of 
the xiab distributions dN / dxi a b for secondary protons produced according to the 
models DPMJET-II.5 and QGSJET. The distributions are plotted for E lab = 
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Fig. 3. (a) Average energy fraction for charged pion ir + , ir~ production K n + n - in p- 
N collisions as function of the (nucleon-nucleon) cms energy <Js. (b) Average energy 
fraction for net baryon B — B production K B _g in p-N collisions as function of the 
(nucleon-nucleon) cms energy y/s. 
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Fig. 4. (a) Spectrum weighted moments for charged pion production f n + n - in p-N 
collisions as function of the (nucleon-nucleon) cms energy \fs. (b) Spectrum weighted 
moments for charged Kaon production f K + K ~ in p-N collisions as function of the 
(nucleon-nucleon) cms energy ^fs. 



10 13 , 10 16 and 10 20 eV . At the lowest energy the two models agree roughly, 
but there are considerable differences in the extrapolations to 10 16 and 10 20 eV. 
The differences between DPMJET-II.5 and SIBYLL-1.7 are mainly due to the 
baryon stopping diagrams implemented in DPMJET-II.5 [Q,^| but not present in 
SIBYLL-1.7. The scaling behavior of the distributions according to DPMJET- 
II.5 and SIBYLL-1.7 is nearly perfect in the xi a b region from 0.1 to 0.6, but this 
is just the region where the non-scaling appears in QGSJET. The behavior of 
the models in p-p and p-p is rather similar to the one in p-N collisions. Clearly, 
measurements of leading baryon distributions at the TEVATRON and at the 
LHC will resolve the conflict between the three models. 

Similarly, in Figs. ||.a and [|b we present the extrapolations to high energy 
Eiab of the xf distributions F (xp) for secondary 7r + produced according to the 
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models DPMJET-II.5 and QGSJET. The distributions are again plotted for E lab 
= 10 13 , 10 16 and 10 20 cV . SIBYLL 1.7 behaves similar to DPMJET-II.5 but 
as seen before DPMJET-II.5 and QGSJET differ considerably in the shape of 
the distributions, only measurements at the highest possible energies can resolve 
the situation. 




Feynmam 



Feynmam 



Fig. 5. (a) a^af) -distributions of secondary protons produced in p-N collisions at Ei a b 
= 10 13 , 10 16 and 10 20 eV according to the DPMJET-II.5 model. (b):r; a f,-distributions 
of secondary protons produced in p-N collisions at Ei a b = 10 13 , 10 16 and 10 20 eV 
according to the QGSJET model. 




Feynmam xp Feynmam xp 

Fig. 6. (a) ^-distributions of secondary 7r + produced in p-N collisions at Ei a b = 
10 13 , 10 16 and 10 20 eV according to the DPMJET-II.5 model. (b)xF-distributions of 
secondary 7r + produced in p-N collisions at Ei a b = 10 13 , 10 16 and 10 20 eV according 
to the QGSJET model. 



In the next four Figures we will demonstrate, that the differences in the 
extrapolations to high energy in Fe-N collisions are at least as large as in p-p 
or p-N collisions. However we find also, that in Fe N collisions QGSJET does 
not run up to such high energies as DPMJET and SIBYLL. 
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For a better comparison with our p-N Figures, we plot also the Fe-N com- 
parisons up to y/s = 1000 TeV per nucleon-nucleon collision. What is needed 
for cosmic ray simulations are lab energies up to 10 20 to 10 21 eV per nucleus, 
the energies reached with QGSJET are sufficient for this. 

In Fig. ^.a we compare average < p± > of produced charged hadrons as 
function of the (nucleon-nucleon) cms energy y/s in Fe-N collisions. The differ- 
ences between the three models look very similar to the differences found in p-N 
collisions in Figjlja. In Fig. |^.b we compare the multiplicity of produced charged 
hadrons as function of the (nucleon-nucleon) cms energy y/s in Fe-N collisions. 
Here the differences between the three models look completely different than in 
p-N collisions in Fig.||.a. 

In Fig. ||,a we compare the average energy fraction for charged pion 7r + , 7t~ 
production w - in Fe-N collisions as function of the (nucleon-nucleon) cms 
energy y/s. At low energies K w + „- rises with y/s for all three models and the 
absolute magnitudes differ by less than 20 percent. However, starting at y/s = 
5 TeV for QGSJET K„+ „- starts to decrease strongly, while for the other two 
models K n + i7r - continues to rise with y/s. 
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Fig. 7. (a) Average transverse momenta for charged hadron production in Fc-N col- 
lisions as function of the (nucleon-nucleon) cms energy y/s. (b) Average multiplicity 
of charged hadron production in Fe-N collisions as function of the (nucleon-nucleon) 
cms energy y/s. 



In Fig. [8|.b we observe a quite similar behavior of the spectrum weighted 
moments for charged pion production f- K + ! - K - in Fe-N collisions as function of 
the (nucleon-nucleon) cms energy y/s. Now the curves for DPMJET-II.5 and 
SIBYLL-1.7 rise in a similar way with y/s , but the absolute size differs consid- 
erably. For QGSJET we find , similarly to the last Figure a decrease starting at 
about y/s = 5 TeV. 
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Fig. 8. (a) Average energy fraction for charged pion 7r + ,7r~ production K w + n - in 
Fe-N collisions as function of the (nucleon-nucleon) cms energy sfs. (b) Spectrum 
weighted moments for charged pion production f 7r + 7T - in Fe-N collisions as function 
of the (nucleon-nucleon) cms energy y/s. 



Comparison of the models after simulating the Cosmic 
Ray cascade 

Next we present one comparisons from the Karlsruhe code comparison J48| , 
which refers to the Cosmic Ray cascade simulated with different event generators. 





Contours in the log 10 iV M -log 10 N e plane for p and Fe induced showers of 
and 10 15 eV . The distributions were calculated using the CORSIKA shower 
with 5 different event generators for the hadronic interactions. The HPDM 
not mentioned before in this talk is a DPM similar model. 
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The distribution chosen in this comparison is motivated by the interest of the 
KASKADE[5fl experiment in Karlsruhe. 

In Fig. | Fe and p induced showers with energies of E = 10 14 and 10 15 eV 
are plotted in the log 10 -log 10 N e plane (Muon-number -Electron-number 
plane). The distribution of events according to each of the 5 interaction models 
for each energy and primary particle is indicated by contours. Considering such 
a plot calculated with only one of the models, where Muon number is plotted 
over electron number, the impression is, that a simultaneous measurement of 
Muon-number and Electron number allows to determine the primary energy as 
well as the composition of the primary component. In such a plot we would 
see, that for instance VENUS and DPMJET agree very well, but the contour 
according to SIBYLL for Fe projectiles of E = 10 15 eV overlaps the VENUS 
and DPMJET contours for p projectiles at this energy. From these differences 
between the models one can conclude, that at present the systematic errors in 
the cascade calculation prevent to identify safely the composition of the primary 
component from such measurements. 



4 Standardization of Monte Carlo hadronic production 
models 

There were in the past and there are at present some activities to standardize 
some features of the Monte Carlo codes. Examples are: 

(i) Universal code numbers for particles and nuclei in the computer codes. The 
Particle Data Group proposed such a Monte Carlo particle numbering scheme, 
see for instance |58| . This scheme is used for the output of events by many event 
generators (examples are PYTHIA, JETSET and other Lund codes, HERWIG 
and the DPM event generators DPMJET and PHOJET), but the scheme is not 
suitable for the internal running of the codes, therefore each event generator 
uses a second (internal) particle numbering scheme. It would certainly be of 
advantage to standardize also these internal particle codes. 

(ii) A standard COMMON block for the presentation of the Monte Carlo 
events was proposed by working groups for the LEP and LHC colliders. The 
structure of this common block has been suggested in Refs. |^9],[50) . This stan- 
dard COMMON follows closely a scheme first used in Monte Carlo codes of the 
Lund group. This COMMON contains not only the final particles produced in 
the collisions, it also documents all the internal features of the model, for in- 
stance: original projectile and target particles, partons, chains formed out of the 
partons, hadronization of chains, decay of unstable particles, intranuclear cas- 
cade interactions, excited nuclei, evaporation particles and stable residual nuclei 
and nuclear fragments. 

(iii) A similar effort has been started by the community of heavy ion experi- 
ments for the event generators for heavy ion collisions |El[] . 
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5 Conclusions 

Most hadronic event generators which can be used for simulating hadronic and 
nuclear collisions up to the highest energies are quite similar in their construc- 
tion and in the underlying theoretical concepts. At energies, where data from 
accelerator and collider experiments are available the models agree rather well 
with each other and with the most important features of the data. As soon as 
we compare the extrapolations of the models at higher energy we find in spite 
of the similarities in the underlying theoretical concepts quite often striking dif- 
ferences between the predictions of the models. We conclude: (i) Measurements 
of inclusive hadron production at the TEVATRON Collider and in the future 
at the LHC collider are very important to guide the models, (ii) More theoret- 
ical efforts to improve the hadronic and nuclear collision models are needed to 
get better extrapolations, (iii) Cosmic ray experiments should use in the high 
energy region simulations with more than one hadronic interaction model for 
interpreting their data. 
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